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ABSTRACT 

How accretion proceeds around the massive black hole in the Galactic center and other highly 
sub-Eddington accretors remains poorly understood. The orbit of the S2 star in the Galactic 
center passes through the accretion disk of the massive black hole and any observational 
signature from such interaction may be used as an accretion probe. Because of its early stellar 
type, S2 is expected to possess a fairly powerful wind. We show here that the ram pressure of 
the accretion disk shocks the stellar wind fairly close to the star. The shocked fluid reaches a 
temperature of ~ 1 keV and cools efficiently through optically thin, thermal bremsstrahlung 
emission. The radiation from the shocked wind peaks around the epoch of the pericenter 
passage of the star at a luminosity potentially comparable to the quiescent emission detected 
from Sgr A*. Detection of shocked wind radiation can constrain the density of the accretion 
disk at a distance of several thousands of gravitational radii from the black hole. 

Key words: accretion, accretion discs — black hole physics — galaxies: active — radiation 
mechanisms: thermal — shock waves — stars: winds, outflows 
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1 INTRODUCTION 

The compact radio source Sgr A* is believed to mark the location 
of the massive black hole in the center of our Galaxy. The Galac- 
tic center is also observed as an IR and, possibly, X-ray source 
of luminosity ~ 10 36 erg s -1 and < 3 x 10 33 er g s -1 , respectively 
Genzel et al.ll2010h. Both IR jGhez et alj|200j) and X-ray flaring 



BaganofF et alj|200lh is regularly observed on timescales ranging 
from minutes to hours. Flares are believed to be associated with 
processes taking place close to the black-hole horizon. 

The radiation observed from Sgr A* is believed to be pow- 
ered by the accretion process. The radiative power is, however, 
a small fraction of the rate at which gravitational energy is re- 
leased by the gas accretion. The accretion is, therefore, expected 
to take place through a hot, quasi-virialized, thick disk forming 
a "Radiatively Inefficient Accretion Flow" or RIAF. The fate of 
the non-radiated energy is model-dependent. The energy may be 
advected into the black ho le through an "Adve ction-Dominated 
Accretion Flow" or ADAF jNaravan et al.|[l995l) ; carried by con- 
vective motions in a "Convection-Dominated Accre tion Flow" or 
CDAF JOuataert & Gruzinovll2000l : iBall et al.ll200ll) : or in kinetic 
form through wi nds from the disk as in the "Inflow-O utflow Solu- 
tions" or ADIOS l lBlandford & Begelmanll999ll2004l) . The various 
models make distinctly different predictions for the gas density and 
its radial profile in the disk, but convincing observational probes 
are still lacking. X-ray observations constrain the electron density 
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close to the sphere of influence of the black hole (i.e., at the Bondi 
radius R b =s 0.04 pc = 2 x 10 s R g , where R g = GM BH /c 2 is the 
gravitational radius for a black hole mass Mbh - 4.3 x 10 s M Q ) to 
be n b ~ lOOcirr 3 jBaganoff et al.ll2003h . but offer little clues for 
the gas properties within that radius, where the accretion disk is 
located. 

The region of R < R\, is not devoid of sources (besides Sgr 
A*). It is filled with tens of massive star s, the so-called S cluster 
(Genzel et al. 20031: rGillessen et"al]|2009h . Most of these stars are 
B dwarfs and have elliptical orbits bringing them as close as ~ a 
fewx 10 3 ^ g from the black hole. B stars are also known to have 
powerful winds of substantial kinetic luminosities L„. > 10 34 erg s 
and characteristic mass loss rates of ~ 10~ 7 M o yr . Among the 
S-cluster stars, the S2 is characterized by both a close pericenter 
passage to Sgr A* and the earliest stellar type, possibly connected 
to the most powerful wind in the cluster. As we show here, the 
interaction of the wind of S2 with the accretion disk is expected 
to result in a characteristic rise of the X-ray emission from Sgr A* 
around the epoch of pericenter passage on a timescale of months, 
that can be used to probe the gas properties at the Galactic center. 



2 THE S2 STAR AND THE ACCRETION DISK IN THE 
GALACTIC CENTER 

The S2 star is th e brightest of the S c luster. It has a ~16 year orbit 
around Sgr A* (Gillessen et al. 2009) of eccentricity e = 0.88 and 



a pericenter distance of R p = 2800 R„. 
the star reaches v p at J2R^/R^ c at 8 



At pericenter, the velocity of 
The nature of 
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the S2 star has been revealed in the study bv lMartins et al.l ( [2008b . 
S2 is an early B dwarf star (of spectral type B0-2.5V) of intrinsic 
luminosity log L/L Q = 4.2 -r 4.8, effective temperature T = (20 -r 
30) x 10 3 K, and radius R, = 10 R . The lack of particular wind 
absorption signatures in its spectrum places an upper limit on the 
wind mass loss rate of M < 3 x 10~ 7 v 8 M Q yr -1 , where v w = 10 8 v g 
cm s is the speed of the wind . 

The iMartins et al. I d2008h upper limit on the wind mass loss 
from S2 is only modestly constraining. Early B stars have, as a 
population, powerful and fast winds. Their wind velocity ranges 
within v„. = (0.5 -=- 2) x 10 s cm s -1 , while their mass loss rate is i n 
the range M = (0.3 -=- 3) x 10 7 M s yr 1 jKudritzki & Prnll2OO0l) . 
Here we adopt v„. = 10 8 cm s and M = 10~ 7 M Q yr -1 as refer- 
ence values for the wind from S2, with the understanding that v„ is 
uncertain by a factor of ~ 2 and M by a factor of ~ 3. 

The accretion rate of Sgr A* is quite uncertain and may well 
be distance-dependent. Theoretical modeling and Faraday rotation 
measurements place limits in the range M BH = 3 x 10~ 9 -r 10~ 7 M 
yr -1 for the accretion rate at the black hole dYuan et al.l 120031 : 
iMarrone et ail |2007| ; iMoscibrodzka et al.l l2009h - Given the esti- 
mated accretion rate, a radiatively efficient disk would radiate away 
L ~ O.IMbhC' 2 > 10 3 7 erg s~'. The bol ometric luminosity from Sgr 
1 erg s -1 



A* is ~ 10 36 erg s -1 jGhez et al.l r2004) suggesting that the most of 



the gravitational energy released is likely not radiated away. Such 
an accretion flow is termed as RIAF. RIAFs are characterized by a 
thick H/R ~ 1 disk that is partially pressure and partially rotation 
supported. The total pressure (sum of the ram and gas pressures) 
of the disk at distance R from the black hole is P ~ pc 2 (R g /R), 
where p(R) is the gas mass density. This expression holds rather 
independently of the details of the adopted model for the flow 
(ADAF, CDAF, or ADIOS). The electron density at the Bondi ra- 
dius R b ~ 2 x l0 5 R g is constrained directly by X-ray observations 
to be n b ~ 100 cnr 3 (Baganoff et al. 2003). The gas density closer 
to the black hole is essentially unconstrained. Possible scalings for 
the density profile adopted h ere are n = n b (R b /R ) 3/2 as, e.g., moti- 
vated by the ADAF solution jNaravan et al J 19951); or n = n b (R b /R) 
as, e.g., motivated by GRMHD sim u lations jlVIcKinnev et alj2012l ; 
iTchekhovskov & McKinnevl |2012| ; iNaravan et alj |2012|) . In the 
case of a CDAF, the den sity profile is shallower n oc R~ l/2 
( Ouataert & Gruzinovl2 000). At the pericenter of S2 (i.e., R = R p ~ 
0.01.Kb), the density of the disk material is, therefore, expected to 
ben- n b (R b /R p ) ll2+3/2 ~ 10 3 -=- 10 5 cm- 3 . 



3 THE STELLAR WIND INTERACTION WITH THE 
ACCRETION DISK 

The interaction of the stellar wind with the accretion disk has two 
potential observational signatures. The stellar wind is terminated by 
a strong shock because of the confining pressure of the disk. The 
shocked gas reaches a temperature of ~ 1 keV and cools via ther- 
mal bremsstrahlung emission mainly in the X-ray band. The RIAF 
also undergoes a shock upon interaction with the stellar wind. The 
shock may accelerate non-thermal particles dNaravan et al.l 12012] ; 
ISadowski et al.l 120131) . Inverse Compton scattering of the bright 
stellar photon field by energetic electrons can power detectable hard 
X-ray and y-ray emission. 



3.1 Thermal emission from the shocked wind 

While the stellar wind expands, its ram pressure drops with dis- 
tance as P v = Mv n /4/rr 2 , where lower case r is measured from 



the center of the star. The stellar wind is terminated by a strong 
shock at a distance r sb where the wind ram pressure is balanced 
by that of the RIAF. The pressure of the RIAF is the sum of the 
thermal pressure and the ram pressure resulting from the relative 
motion of the RIAF and the star. The thermal pressure is of order 
Pth ~ pc 2 (^ g /2^) while the ram pressure is mainly a result of the 
stellar motion, so that P ram ~ 2pc 2 (R e /R). The ram pressure varies 
by ~ 50% depending on the angle of the stellar and disk orbits. 
For the purpose of our estimates, we set the total disk pressure as 
P tot ~ 2.5pc 2 (R s /R). Equating P tot = P w at R - R p gives the dis- 
tance at which the stellar wind is shocked, when S2 is at pericenter: 



r sh = 6 x W l3 M l ' 2 vl' 2 n: 112 cm, 



(1) 



where M = 10~ 7 M_7 M Q yr~' is the stellar wind mass loss, and 
n = 10 4 «4 cirr 3 is the electron (or proton) number density of the 
RIAF at pericenter. Note that for reasonable parameters r sh <sc R p , 
so the wind is terminated at a "small" distance from the S2 star. 
In the post-shock region, the temperature of the gas is T sb = 

1.3 x 10 7 v^ K, and the electron number density is n sb = 4n w (r sb ) = 
3 x 10 6 M4Vg 2 cnr 3 , where n n . is the electron number density in the 
unshocked wind. Here, we have assumed solar metallicity (our best 
guess for the wind from a main sequence star). The wind cools 
through optically thin, thermal bremsstrahlung emission. For tem- 
peratures T ~ 10 7 K and a solar-like composition, line emission 
dominates the cooling rate over free-free emission by a modest fac- 
tor. The emissivity is approximately constant for T ~ (0.3 -r 3) x 
10 7 K and of order A N ~ 3 x 10~ 23 erg cmV 1 dStevens et al.ll 19921; 
ISutherland & Dopitall993 j ). The cooling timescale of the plasma is 
then f c = 2kBT sb /n sb A^ = 4 x lCPv^n^' s. bremsstrahlung cool- 
ing competes with the adiabatic expansion of the shocked wind. 
The latter takes place on a timescale f exp ~ r s hA'sh ~ 4r sh /v w — 

2.4 x 10 6 M^ 2 Vg 1/2 n^ 1/2 s, where v sh = v w /4 is the post-shock flow 
velocity^] Out of the total kinetic luminosity L„. = Mv 2 ,/2 of the 
wind, a fraction f cxp /? c - 0.06 MjjnV 2 vZ 9 is radiated in the X- 
ray band as thermal bremsstrahlung. The X-ray luminosity of the 
shocked wind region is 



— L„. 



2 x 10" M 



33^3/2 n 1/2 v -5/2 



(2) 



Given that the quiescent X- ray luminosity from S gr A* is at a sim- 
ilar level (a few 10 33 erg/s; Ba ganoff et ail 2003). it is possible to 
measure contributions to the total emission from the shocked wind 
when the S2 star is close to pericenter. 

In Fig.[T] we present the X-ray light curves expected from the 
shocked wind of S2 when the star approaches its pericenter, in mid- 
2018. For each point along the stellar orbit - characterized by the 
inclination i and the argument of periapsis a> - we extract the lo- 
cal conditions in the disk from the model bv lSadowski et al. I d2013l 
see Table Al there), which is based on GRMHD simulations of the 
accretion flow around Sgr A*0 We compute the shock radius r sh at 
each time by balancing the wind ram pressure with the total (ram 
plus thermal) pressure in the accretion flow. When calculating the 



The expansion time of the shocked wind f cxp is somewhat shorter than 
the typical evolution time ? cv = R p /v p of the stellar orbit at pericenter (i.e., 
the time over which the confining disk pressure changes). We, therefore, 
assume that the wind-disk interface evolves quasi-steadily along the stellar 
orbit. 

2 The inclination i is such that ;' = 0° if the orbital plane coincides with the 
disk equatorial plane. The argument of periapsis to is measured along the 
orbital plane from the line of nodes, where the orbital plane intersects the 
disk midplane. We have ui = 0° if the pericenter is on the disk midplane. 
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Figure 1. Lightcurves of the X-ray emission from the shocked stellar wind, for different stellar and disk parameters. In all the panels, the vertical black line 
marks the time of pericenter. (a) We vary the inclination of the stellar orbit with respect to that of the accretion disk (i = 0° in red, (' = 30° in green, i = 60° in 
blue), for a fixed argument of periapsis (to = 90°). Solid lines show counter-rotating orbits, dashed for co- rotating orbits, (b) We vary the argument of periapsis 
of the stellar orbit (at = 0° in red, to = 45° in green, to = 90° in blue), for a fixed orbital inclination (i = 60°). We only show counter-rotating orbits, (c) We 
vary the density profile in the accretion flow (n oc R^ in red, n oc Rr^ 1 in green, n oc R~ 3 ' 2 in blue), for a fixed density at the Bondi radius tit = 130 cirT 3 . We 
only show counter-rotating orbits with (' = 0°. (d) We vary the properties of the stellar wind with respect to our standard choice M-j) = (1,1) plotted in 
red, by changing the wind velocity (green) or the mass loss rate (blue). We only show counter-rotating orbits with i = 0°. 



disk ram pressure as the star moves through the disk, we differen- 
tiate between stellar orbits that are co-rotating vs counter-rotating 
relative to the disk azimuthal velocity. Given the shock radius r sb at 
each time, we compute the X-ray luminosity as outlined above. 

In panels (a) and (b), we show how the X-ray emission de- 
pends on the orientation of the stellar orbit relative to the disk plane 
(we vary the inclination i in panel (a), for a fixed argument of pe- 
riapsis to = 90°; and we change the argument of periapsis to in 
panel (b), for a fixed inclination i = 60°). As shown in panel (a), 
counter-rotating orbits (solid lines) always result in stronger X-ray 
signatures than their co-rotating counterparts (dashed lines, with 
the same color), since the ram pressure from the disk confines the 
shock closer to the star, giving a higher post-shock particle density 
and enhanced bremsstrahlung emission. The difference in the peak 
luminosity is, however, modest (at most a factor of 2) when com- 
paring counter and co-rotating orbits. The X-ray emission peaks 
when the stellar orbit intersects the disk midplane, since the disk 
density is the highest there (and so the shock is closest to the star). 
Depending on the argument of periapsis, the star may cross the 
disk only once, at pericenter (red line in (b), for to = 0°); or two 
times, symmetric around the pericenter epoch (blue line in (b), for 
to = 90°); or two times, but asymmetric with respect to pericenter 



(green line in (b), for to = 45°). The shape of the light curve around 
pericenter can then be used to constrain the orientation of the orbit 
of S2 relative to the accretion flow of Sgr A*. 



The peak luminosity depends significantly on the struc- 
ture of the accretion flow around Sgr A* and on the prop- 
erties of the wind of S2. In panel (c), we fix the disk den- 
sity at the Bondi radius n b = 130 cm -3 and we compare the 
scaling n oc R ^ exp ected on the basis of GRMHD simula- 
20121; iTchekhovskov & McKinne\l 12012); 



tions (McKinnev 



R exp ' 
ey ct al. 



iNaravan et al .1120121 ; Sadowski et al.ll2013l) with the prediction n i 
Rr i!1 of the ADAF solution (in blue) and with the profile n oc R~ 1/2 
of CDAF models (in green). In panel (d), we vary the velocity and 
mass loss rate of the stellar wind, with respect to our reference 
choice (vg,M_7) = (1, 1), plotted as a red curve. The trends seen 
in panels (c) and (d) can be simply understood from the scalings il- 
lustrated in eq.[2] which gives the luminosity expected at pericenter. 
Once the properties of the stellar wind of S2 are better determined, 
the X-ray emission from the shocked wind could be used to place 
important constraints on the density profile of the accretion flow 
around Sgr A*. 
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3.2 Non-thermal radiation from the shocked RIAF 

ISadowski et al, recently showed that a highly elliptical orbit 

of an object moving through a RIAF is likely to be supersonic for 
the disk itself. The Mach number of the shock is 1.5 -r 3 depending 
of the relative orientation of the stellar and disk plane. The RIAF is 
shocked and elect ron acceleration has b een proposed to take place 
at the shock front dSadowski et al.ll2013f) . We explore this possibil- 
ity by assuming that a fraction e ~ 0.1 of the total dissipated energy 
is deposited into non-thermal electrons with a rather flat power-law 
energy distribution (i.e., an energy spectrum dN/dy oc y~ 2 ). Rel- 
ativistic electrons will produce X-rays and y-rays by upscattering 
the powerful UV radiation from the S2 star. 

The RIAF material approaches the star with speed v ~ v p , dis- 
sipating kinetic energy at a rate La ss ~ 7ir 2 h pv 3 ~ 10 35 M_7Vserg s _1 . 
The shocked RIAF material expands on a timescale f cxp ~ 4r sh /v p a 
3 x 10 5 M [ J 2 Vg 2 n 4 ' 12 s. A relativistic electron with Lorentz factor y 
cools on a timescale t c = 3 x 10 7 /yl/ p h s, where l/ p h = L,/4nr^ h c is 
the energy density of the stellar radiation field, and we assume that 
the S2 star has luminosity L, = 10 38 L» 38 erg s" 1 . The characteris- 
tic energy of the stellar photons is E, ~ 3k^ T, ^ 6 eV, requiring 
electrons with Lorentz factor y = 70 E 1 ^^ for the production of ~ 
30keV photons. A fraction t cxp /t c = 0.05 L.^E^MZ^v^' 2 /! 1 / 2 
of the total dissipated power Ldiss is radiated away. The luminosity 
of the IC component at photon energy E is 

L,c ~ e^L diss ~ 6 x \Q^e^L, M E^ keV M[ l2 vfnf erg s" 1 . (3) 

Currently the quiescent level of emission in the hard X-rays from 
Sgr A* is not known. Extrapolation of the soft X-ray emission spec- 
trum indicates that the level of hard X-ray luminosity at ~30keV 
may be in the range ~ 10 32 -r3x 10 33 erg/s. We conclude that L ic may 
be powerful enough to be detectable in this band over the quiescent 
emission from Sgr A*. 

Electrons with y > 10 3 M U2 v 1 ^ 2 L^n 4 i/2 are in the fast cool- 
ing regime (r cxp /f coo i > !)■ As a result, all the energy that is in- 
jected in those electrons is radiated away in the form of E > 
6 M-jVgL^n^MeV photons. Under the optimistic assumption of 
a flat injection electron distributior0 up to y > 10 4 , as much as 
Lie ~ 10 34 e_iM_7Vgerg/s may be radiated in the MeV band [and 
up to ~ 200(6eV/£,)GeV, where Klein-Nishina suppression ap- 
pears]. Unfortunately, the high density of sources in the Galactic 
center region makes detecting such a y-ray signal very challenging. 



4 DISCUSSION AND CONCLUSIONS 

Given its stellar type, S2 is expected to posses a moderately pow- 
erful stellar wind with M ~ 10~ 7 Mq yr~' and v, v ~ 10 8 cm s -1 , 
Upon interaction with the accretion flow of Sgr A*, the stellar 
wind is shocked fairly close to the star, resulting in substantial 
bremsstrahlung cooling in the ~a few keV band. The predicted 
emission spectrum is that of optically thin, thermal bremsstrahlung, 
characterized by strong emission lines. The peak of the emission is 
predicted to take place around the pericenter passage of the star 
and the emission should remain bright over a several-month period 
(Fig. 1). The X-ray luminosity (eq. 2) is potentially detectable close 
to its peak (or peaks) as an increament of the quiescent emission 

3 If the electron power-law slope is p > 2, the IC luminosity at a photon 
energy E will be suppressed by a factor of ~ (E/E,)( p ~ 2 ^ 2 with respect to 
the estimates given above. 



that lasts for monthfl The slow evolution pattern of the emission 
from the wind-disk interaction should allow to discriminate it from 
the observed rapid (minute to hours) flaring that is probably related 
to activity close to the black-hole horizon. The shape of the light 
curve around the peak(s) depends on the inclination or the orbit of 
S2 with respect to the disk and can be used to determine the orbital 
plane of the disk. 

The wind-disk interaction leads to the formation of a shock in 
the disk material as well. If the shock of the RIAF accelerates non- 
thermal electrons, then IC scattering of the stellar radiation field 
leads to a significant hard X-ray signal. NuSTAR dHarrison et al.l 
l2013h has the sensitivity to detect such a hard re-brightening from 
Sgr A* at the epoch of the pericenter passage of S2. 

Similar considerations apply to the wind-disk interactions of 
other S stars. From eq. 2, it is clear that high wind mass loss rates, 
slow wind velocities and small pericenter radii result in the bright- 
est events. However, the other members of the S cluster are dimmer 
than S2 and are likely to have weaker winds. For the moment S2 
remains the best candidate for such a study. 

The next pericenter passage of S2 will take place in 2018. By 
that time, a better observational understanding of the quiescent lu- 
minosity of Sgr A* may be in place in both soft and hard X-rays. 
This will facilitate the detection of even weak enhancements of 
the X-ray luminosity of Sgr A* on a timescale of months to years 
around the closest approach of S2. Furthermore, our work moti- 
vates the need for a more precise determination of the wind prop- 
erties of S2 (such as mass loss rate, velocity and metallicity). The 
curre nt limit for the wind mass loss rate from S2 is M < 3 x 10~ 7 M o 
yr _1 dMartins et alj|2008h . which is only modestly constraining. A 
precise assessment of the wind properties will make any detection 
of shocked wind emission (or even upper limits) a very powerful 
probe of the accretion disk. 

In the Summer of 2013, the G2 cloud dGillessen et~aUl2012l 
l2013l) is passing at its pericenter in the disk of Sgr A* at a distance 
of ~ 5000 R g . Interactions of the c loud with t he accretion disk may 
lead to X-rav dGillessen et al.ll2012l) or radio dNaravan et alj|2012l : 
Sado wski et al.ll2013h signatures that can be used to probe the disk 
properties. An interesting fact is that the orbit of the S2 star is 
highly inclined with respect to that of the G2 cloud, providing a 
potential probe of the disk properties on a different plane. 
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